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Abstract. Lie's linearizability criteria for scalar second-order ordinary differential equa- 
tions had been extended to systems of second-order ordinary differential equations by using 
geometric methods. These methods not only yield the linearizing transformations but also 
the solutions of the nonlinear equations. Here complex methods for a scalar ordinary dif- 
ferential equation are used for linearizing systems of two second-order ordinary and partial 
differential equations, which can use the power of the geometric method for writing the solu- 
tions. Illustrative examples of mechanical systems including the Lane-Emden type equations 
which have roots in the study of stellar structures are presented and discussed. 



1. Introduction 

Lie developed a systematic procedure for solving nonlinear ordinary differential equations 
(DDEs) with some minimal symmetry under transformations of the dependent and indepen- 
dent variables, called point transformations, by using group theory [22] . For second-order 
ODEs he provided criteria for their being transformable to linear DDEs, provided they are 
maximally symmetric [23]. The requirement is that they be at most cubically nonlinear 
in the first derivative and that their coefficients satisfy a set of four constraints for the 
first derivatives, involving two auxiliary functions. Tresse [371 EH] reduced the number of 
constraining equations to two for higher derivatives by eliminating the auxiliary functions. 

Since Lie's time there have been various developments following up his work. There was 
the development of contact transformations (see e.g., [TJ [TS]) and their use for linearizing 
two classes of third-order ODEs [TU [17] ; approach of Cartan type used for the same purpose 
[T0| [TT] : a method for a special case of a class of third-order scalar ODEs [32]; Lie's alge- 
braic methods used for classification of systems of ODEs p5l |26l [39]; extension to partial 
differential equations (PDEs) using potential symmetries and otherwise [3 [Gj [21]; the use 
of algebraic computing to extend Lie's methods to general third- and fourth-order ODEs 
[T^ I2U| [23]; conditional linearizability [2H1 [21]; geometric methods for symmetry analysis 
[U [131 [3S] and their use for linearization [301 [31] ; and the development of complex symmetry 
analysis [H [21 [3]. Here we use the last two developments to provide linearization for a class 
of systems of two PDEs and an alternative method for linearizing a class of systems of two 
ODEs that is not equivalent to the earlier method for the same purpose [3J. 

The geometric approach for linearization of a system of two ODEs requires that the system 
be at most cubic in its first derivatives and satisfy a generalized set of Lie- Tresse invariant 
conditions which are written in terms of the coefficients of the system of equations [301 131] • 
This requirement comes from regarding the system as a projection of the geodesic equations 
in a fiat space (in curvilinear coordinates). The linearizing transformation is then obtained 
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by converting the metric in the given coordinates to one in Cartesian coordinates. This en- 
ables us to write down the solution of the nonlinear equations directly. In complex symmetry 
analysis one studies the relationship between the algebraic properties of the complex differ- 
ential system and the corresponding real system of ordinary or partial differential equations, 
including the Cauchy-Riemann equations, which arises from the complex system via com- 
plex splitting of the dependent and independent variables PQEl E] . We present linearization 
criteria, both invariant and algebraic, for systems of two real ODEs and PDEs that arise 
from the linearization of scalar complex second-order ODEs. This also provides an answer 
to the inverse problem of those systems of ODEs and PDEs that arise from the linearization 
of scalar complex second-order ODEs. Furthermore, we provide examples of mechanical 
systems that are linearizable by our procedure. 

The plan of the paper is as follows. In the next Section the salient points of complex 
symmetry analysis and the geometric methods developed are mentioned. In Section 3 invari- 
ant linearizability criteria for systems of two PDEs are given and some physical examples 
provided. In the subsequent section restricted complex transformations is introduced and 
used to yield invariant linearizability criteria for systems of two ODEs. Again, illustra- 
tive examples of mechanical systems are provided. Section 5 consists of a discussion and 
conclusion. 

2. Preliminaries 

A scalar ODE for a complex function of a complex variable can be written in terms of the 
real and imaginary parts of the dependent variable as functions of the real and imaginary 
parts of the independent variable. As such it yields a pair of PDEs. Of course we need to 
include the Cauchy-Riemann equations (C-REs), which guarantee the differentiability of the 
dependent complex variable in our system. As such we get a system of four PDEs, two of 
which are first-order equations. The symmetries of the original ODE are not identical with 
those of the system of PDEs. Nevertheless the solutions of the ODE obtained by using its 
symmetries give the solution of the PDE. As stated this may seem trivial. However, one 
could pick a system of PDEs and check if it corresponds to an ODE. In the present note we 
also specify those systems of PDEs that could be transformed to such complex differential 
equations. Then one can use the solution of the ODE to write a solution of the system of 
PDEs. The correspondence between the symmetries of the system of PDEs and the ODE is 
the subject of "complex symmetry analysis" [Ill2]- Thus, for example, the single infinitesimal 
generator for a function of one variable, = d/dz becomes the pair of generators X. = d/dx 
and Y = d/dy, where z is the complex independent variable, x is its real part and y its 
imaginary part. Again the scaling symmetry for the independent complex variable gives the 
scaling symmetry for the two real variables together along with the rotation generator in 
the 2-dimensional space. Similar remarks apply for the generators involving the dependent 
complex variable or a mixture of the two. 

For a linear scalar ODE the conversion to the complex form and thence to the system of 
two PDEs is trivial. The only complication is that here the system becomes a set of four 
equations, the two C-REs and the two PDEs for the two real variables corresponding to 
the single complex variable. This can change the symmetry structure but does not make 
any other substantive difference. For a nonlinear ODE the situation changes. Now there 
can be non trivial mixing between the real and imaginary parts of the complex dependent 
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variable so that the system of two PDEs becomes significantly coupled. The question can 
arise whether the C-REs continue to hold under some transformation of the independent and 
dependent variables. This is of relevance for us as we need to use linearizing transformations. 
Suppose the ODE is written for a complex analytic function of a single variable, u{z). Since 
the linearizing point transformation C : {z,u) — )■ {Z,U) is analytic (by definition), the real 
linearizing transformation TZC : (x, y, /, g) — t- (X, Y, F, G) satisfies the C-REs. Note that 
when we go to deal with systems of ODEs, the transformed variables are not guaranteed to 
satisfy the C-REs rather they satisfy a partial analytic structure. In that case we have to 
use the C-REs in the original variables. This creates complications in the linearization of 
those ODEs. Apart from this complexity the solution of nonlinear system is still attained 
through a procedure analogous to analytic continuation. 

The system of geodesic equations is of second-order and quadratically semi linear in the 
first derivatives and has no other terms in it. It inherits the isometries, but can have many 
other symmetries. It was noted |1] that projection of this system, using the translational 
invariance symmetry of the geodesic parameter yields a system of cubically semi linear ODEs. 
The linearizability of a quadratically semi linear system of ODEs of geodesic type is provided 
by regarding the coefficients of the quadratic terms as Christoffel symbols and verifying 
whether the resulting Riemann tensor is zero or not [30j. If the system is linearizable, one 
can find the linearizing transformation by taking the coordinate transformation from the 
metric tensor constructed from the Christoffel symbols p!^ to the metric tensor in Cartesian 
coordinates. The same procedure can be extended to the projected system of cubically semi 
linear ODEs [31]. The power of this method is apparent from the fact that one can not 
only find the linearizing transformations but also be able to write down the solution of the 
nonlinear equation. 



3. Invariant Linearizability Criteria for Systems of PDEs 

First-order scalar ODEs are always linearizable. We only need to deal with second-order 
systems for checking linearizability. The geometric linearization of systems of two ODEs [30| 
|3T] gives a requirement that the system be (at most) cubically semi linear in the dependent 
variables and satisfy a generalization of the Lie conditions coming from the fiatness of the 
space in which the geodesies lie. However, this procedure does not cover all the classes 
of linearizable systems as it requires a 15— dimensional symmetry algebra, whereas there 
can also be 5— ,6— ,7— and 8— dimensional algebras of linearizable systems [2S]. Since the 
Lie symmetry algebra for linear PDEs is infinite dimensional in general, the question arises 
whether we can expect the PDEs to be (at most) cubically semi linear as well. The geometric 
approach relied on the connection between geometry and systems of second-order ODEs via 
the system of geodesic equations [U [13] . As there is no apparent geometrical way of extending 
the Lie conditions to systems of PDEs, it is not clear what the analogues of the linearizability 
criteria would be. 

We extend Lie's linearizability criteria to a class of systems of two PDEs obtainable from 
complex scalar ODEs. The real transformations for linearization of this system of nonlinear 
PDEs can be obtained by decomposing the complex transformations that linearize the com- 
plex ODE. We consider a simple equation to illustrate this. For example, we exploit the real 
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transformation, 



F 



f 



^-x, G 



(1) 



P + 9' 



P + 9' 



to map a nonlinear coupled first-order system of PDEs, 



9x- fy^ -4=f9, 



(2) 



into a linear system 




0, 
0. 



(3) 



Note that the real transformation (1) is equivalent to a complex transformation 



U^{l/u)-z. 



(4) 



Therefore we can employ complex variables to map systems of nonlinear PDEs into their 
simpler forms. It is important to mention that the linearization above is different from the 
classical way of linearizing equations because the real transformation (1) is obtained from the 
complex Lie transformation (4) that may or may not be a real Lie transformation. Therefore 
complex linearizability works in a different way. A natural question following the above 
discussion would be "which systems of PDEs can be treated via complex linearization?" . We 
provide a class of systems of second-order PDEs which upon satisfying a set of four conditions 
is subject to complex linearizability in Theorem 1. It further provides a partial answer to 
another important question, i.e. "which systems of PDEs correspond to complex differential 
equations?" It is indispensable to use complex linearization in several cases because it plays 
a crucial role in extracting solutions of systems of differential equations that would have 
been difficult otherwise. We highlight this feature of complex variables in the examples. In 
the next Section we extend the above characteristic of complex transformations in dealing 
with systems of ODEs. Namely it is shown that, if the above transformation is restricted on 
a single line, then a nonlinear system of first-order ODEs can also be linearized. In general 
restricted complex transformations can be brought into play in the linearization of systems 



For the above purpose we first prove a few results that are used to obtain linearizability 
of this class of systems of PDEs. 

Theorem 1. The class of systems of two second- order partial differential equations 



of ODEs. 



9xx 



XX 



fyy + '^9xy = ^wi{x, y, /, g, h, I), 
9yy - '^fxy = 4w2(x, y, f, g, h, I), 



(5) 



where 



'^h^ f^ + gy, 21 ^g^- fy, 
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(6) 



is complex-linearizable if and only if the functions, wi and W2, are at most cubic in h and I, 
i.e. 



fxx - fyy + '2gxy = 4A^/i^ - l2A^hf - 12A^hH + AAH^ + ABW - ABH^ - 

SB^hl + AC^h - AC^l + AD\ 
Qxx - Qyy - Vxy = I'^AWI - AAH^ + AA^h^ - UA^hf + SB^hl + AB^h"^ - 

ABH'^ + AC'^h + ACH + AD'^ (7) 

together with the constraints on the coefficients 

"iAl^ - ^Aly + QAly + QC^Al + 6CMJ - QAlC^ + 6CMJ - 
QA)D^ - QD^A] + QD^A] - QD^A] + QA^Cl + QA^C^ - 6A^C^ + 
6A^Cl + C]f - + 2C}g - UA^Dj - UA^D] + 12^2^^ - UA^D] + 
2B^C) + 2B^Cl - 2B^C^f + 2B^Cl - AB^Bl - AB^Bl + AB'^BI - AB^B}, - 

^Blf - 2Bl, - 2B% + 2Bl = 0, 
3^L - - ^Al^ + QC'Al + 6C^AJ + GA^C^ - GCUJ - 
^D'^A) - 6D^Al - 6D^Aj + GD^A] + 6A'^Cl + dA^C^ + 6A^C^ - 
6A^Cl + Cjf - - 2C}g - 12A^D) - UA^D] - 12A^D) + 12A^D] + 
2B'^C} + 2B'^Cl + 2B^C} - 2B^Cl - AB'^BI - AB'^B^ - AB^Bl + 

AB'Bl - 2Blf + 2Blf + 2S^^ - 2Bl = 0, 

12D^Ai + 12D^Al - 12D'^Al + UD'^A}, - QD^Bl - QD^Bl + 

y ^ y J y 

eD^Bj - QD'^B] + QA^Dl + QA'-Dl - QA^Dl + GA'DJ + 

Bl^ - + 2B% - 2Cl^ 2Cl^ - 2C% + 20^ - QB'D) - 
QB^D] + QB'^D) - QB^Dl + 3D]j - 3D]g + QD)^ + AC^C] + 
AC^Cl - AC^C} + AC^C] - 2C^Bl - 2C^Bj + 2C^Bl - 20^3] = 0, 
12D^Al + 12D^Al + 12D^Al - UD^A] - QD'^B) - QD'^B^ - 
6D^Bj + QD^B] + QA^Dl + QA^Dl + GA^D^ _ g^i^i ^ 
- C - 25^. - 2C,^ + 2Clf + 2Ci^ + 2C^2^ - QB^D) - 
QB^Dl - GB^Dj + QB^Dl + 3DJ^ - 30^^ - 6D]g + AC^C] + 
AC^Cl + 4C^CJ - AC^C] - 2C^Bl - 2C^Bl - 2C^Bl + 2C^Bl = 0, (8) 

where all the coefficients A\ B\ and D\ (i = 1, 2), are functions of x, y, f and g. 
Proof. We firstly assume an analytic structure on the manifold and that there exists a 
complex transformation that project the functions f{x,y) and g{x,y) to a single complex 
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function u of complex variable z. Furthermore assume that there exists four complex func- 
tions, A, B, C and such that 

A{x,u) = A^{x,yJ,g) + iA^{x,yJ,g), 
B{x, u) = B^{x, y, /, g) + iB'^{x, y, /, g), 
C{x, u) = C\x, y, /, g) + iC'^{x, y, /, g), 

Dix, u) = D\x, y, f, g) + iD\x, y, f, g). (9) 

By invoking (9) we can map the system (7) to a second-order complex differential equation 

u"{x) = A{x, u)u'^ + B{x, u)u'^ + C{x, u)u' + D{x, u), (10) 

which is at most cubic in its first derivative. Therefore it satisfies Lie's linearizability criteria. 
Moreover the conditions (8) can be projected down to a set of two equations 

3A,, + 3A,C + 3AC, - 3A^D + C^u - 6AD^ + BCu - 2BB, - 25,„ = 0, 

QA,D - 3B^D + 3AD, + B,, - 2C,u - 3BD., + 3D„„ + 2CCu - CB, = 0; (11) 

which may be recognized as the Lie compatibility conditions. As the ODE (10) is linearizable, 
so the system of PDEs (7) is also linearizable. 

Theorem 2. // the system of PDEs (5) admits four real symmetries Xi, Yi, X2 and Y2, 
such that 

Xi = P1X2 - P2Y2, Yi = P1Y2 + P2X2, (12) 
for nonconstant pi and p2, and their commutators satisfy 

[Xi,X2] - [Yi, Y2] = 0, [Xi, Y2] + [Yi,X2] = 0, (13) 

then, there exists a point transformation {x,y,f,g) — > {X,Y, F,G), which reduces Xi, 
Yi,X2 andY2 to their canonical form 

d d d d d d 

and the system (5) can be reduced to the linear form 

Fxx - Fyy + 2GxY = 41^1 (X, F), 

Gxx - Gyy - 2FxY = 41^2(X, Y). (15) 

Proof. Suppose that X^ + iYa = Zq, for a = 1,2. Then equation (13) can be replaced 
by [Zi, Z2] = 0, which implies that the two complex symmetries Zi, Z2, commute with 
each other. Further, setting Zi = p{z,u)Z2 for a nonconstant complex function p, justifies 
equation (12). It is proved by Lie (see |21]) that every scalar second-order ODE (10) that 
admits two commuting symmetries such that Zi = p{z,u)Z2 can be transformed into a 
linear ODE U" = W{Q, by applying the point transformation, ( = ({z,u),U = U{z,u) 
which reduces Zi and Z2 to their canonical forms 

d d 

Zi = Z2 = Ctt^- (16) 

The point transformation, {x,y,f,g) — > {X,Y, F,G), can be obtained for complex trans- 
formation which then can be used to convert system (5) into the linear form (15). 
Examples 

1. Higher- dimensional Coupled System of Modified Lane-Emden Type: 
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The Lane-Emden equation arises in the study of stellar structures [9]. We first investigate 
the linearizability of a system of two cubically semi linear PDEs 

fx. - fyy + 2g,y = -12fh + 12gl - Af + 12 fg\ = gy, 
g.. - gyy - 2f^y = -Ugh - 12// - 12 fg + V, /, = -g^. (17) 

The above system can be regarded as a special case of a higher- dimensional Lane-Emden 
system equipped with an analytic structure. The coefficients satisfy linearizability conditions 
(8). Therefore the above system is linearizable. Thus one can check the linearizability 
of a class of systems of PDEs analogously to Lie's technique. In order to construct the 
transformation we employ the symmetries of (17). It admits Lie symmetries 

X - ^ Y - ^ 

OX ay 

^ d d „ d d d d d „ d 

^2 = + 2/^ - - Y2 = y— -x— + /— , (18) 

ox oy of Og Ox Oy Oj Og 

which satisfy all the conditions of Theorem 2. Therefore (18) leads to the linearizing trans- 
formation, 

X = x-—^—r, Y = y+ ^ 



F = -{x-y) — — , G = xy-— , (19) 

2 xf +yg yj - xg 

and we deduce the following system of linear PDEs 

Fxx — Fyy + 2Gxy = 0, 

Gxx - Gyy - 2Fxy = 0. (20) 
2. Consider the nonlinear anharmonic oscillator system given by the PDEs 

fifxx - fyy + 2g^y) - g{g^^ - gyy - 2f^y) = A{h^ - f) - A{f - g'^)wi + 8/5(^2, 
f{.gxx - gyy - 2py) + gifxx - fyy + 2g^y) = 8hl - Sfgwi - 4(/^ - g'^)w2, 

fx gyi fy gxj 

(21) 

where both wi and W2 are arbitrary functions of x and y. We can transform the above system 
into a system of the form (7). Consequently it can be verified that the coefficients of h and 
I satisfy the conditions (8). Therefore the above system is linearizable. Notice that we can 
convert the system (21) into the second-order complex nonlinear anharmonic oscillator ODE 

uu" + w{z)u^ = u''^. (22) 

It admits the symmetries 

Zi=zu—, ^2=u—, (23) 
which yields the complex transformation 

Z = -, U = -\ogu. (24) 

z z 
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The above transformation gives the real transformation 

X —V 

X=— Y ^ 



2-2 _|_ y2 ' _|_ y2 ' 

x'^ + 



F = {x In {f + g') + 2y arctan (g/f)) , 



G = {2x arctan (^//) - y In (/^ + ) , 

X -|- (/ 

that transforms system (21) into the hnear form 

Fxx - Fyy + 2Gxy = -^2^ ((^' " 3^1^')«;i - (Y^ - ^X'Y)w2) , 
Gxx - Gyy - 2FxY = ((^' - 3Xr2)^2 + {Y^ - SX^Y)wi) , 

where 

X -Y 

Wi — Wi 



W2 = W2 
3. We use the transformation 



X -Y 



X2 + y2 ' ^2 + y2 



X = 2/-x2 + |/2, r = 2<7-2a;|/, F = x, G = |/, 
to hnearize the system of PDEs 

fxx - fyy + 2^^j, = 4 + 4 {{h - xf -{I- yf) wi - 8{h -x){l- y)w2, 

Qxx - Qyy - '^fxy = 8(/l - x){l - y)wi + 4 ((/l - xf - {I- yf) W2, 

fx — Qy-i fy ~ ~9xi 

where Wi and W2 are arbitrary functions of the mixed variables 2/ — + y^ and 2g 
The above system possesses the symmetries 

^ d d d d d d 

ox a J og og oj oy 

„ d d . o 2\ 9 ^ d 
^2^x— + y— + {x -y )— + 2xy—, 
ox oy OJ og 

The linearized system of PDEs is 

Fxx - Fyy + 2Gxy = -2{Hwi - LW2), 
Gxx - Gyy - ^Fxy = -2{Hw2 + Lwi). 

These conclude our examples. In the next Section we present another important use 
complex method. 



4. Restricted Complex Transformations and Invariant 
linearizability criteria for systems of two odes 



In the previous section we utilized complex functions of complex variables to obtain lin- 
earizing transformations for systems of PDEs. If we restrict our complex functions to depend 
upon a single real variable, then it would generate different transformations. They trans- 
form systems of ODEs into other systems of ODEs. We extend our three-dimensional space 
of two dependent and one independent variables to a two-complex-dimensional space. In 
the intermediate steps we move off the real line to obtain our linearizing transformations. 
Then the solutions of systems of ODEs are recovered by restricting the independent variable 
to the real line. The procedure is reminiscent of analytic continuation. In the process we 
loose the C-REs so our system is not the same as that with which we started. It is very 
interesting to see the invariance of analytic structure under the restricted complex transfor- 
mation. In this case the ODE is written for an analytic function of a single real variable, 
u{x). It is important to understand how the complex transformation works in the restricted 
domain. The linearizing complex point transformation C : (x, u) — )■ (x, U) can exhibit dual 
nature in terms of its analyticity which yields the following real linearizing transformation 
TZC : (x, /, g) — >• (x, T, (). Because in the transformed variables x can be either complex or 
real. If it is real, then the transformed variables satisfy C-REs only in / and g in which case 
there is a partial analytic structure on the transformed variables, but, if x is complex, then 
the complete analytic structure is restored on the transformed variables. We illustrate this 
important feature in the two examples below. Thus we cannot guarantee that at the end 
a similar linearized system can be obtained via the linearizing transformations from other 
approaches even though we do get the solution by linearization (in the complex). 

The invariance properties of a system of ODEs, 

r = w,ix,f,gj',g'), 

g" = W2{x,f,gJ',g'), (32) 

have been investigated by using complex functions in [H [2] . The idea is to make use of the 
transformation, 

u{x) = f\x) + ig{x), w{x, u) = Wi{x, /, g) + iw2ix, /, g), (33) 

to convert the system (32) into the single ODE, 

u"{x) = w{x,u,u'), (34) 

and then use the standard Lie procedure of linearization. We call this procedure complex 
linearization even though we have no guarantee that the transformed system can be linearized 
via other approaches. To comprehend the connection between complex transformations and 
linearization we firstly revisit an earlier example from the previous section. 
Consider a first-order two-dimensional Riccati system 

f = -f + 9', 

9 = -Vg. (35) 

which is a special case of a general Riccati system in two dimensions. A natural question 
arises: which transformation can linearize the above system? If such a transformation exists, 
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then how can we find it. The use of the transformation, 

r +9 f + 9 

maps system (35) into the simplest system 

r = 0, C' = 0. (37) 

The transformation (36) is a mere consequence of the same complex transformation (4) that 
we used to hnearize system of PDEs (2) in the remaining Section. The only difference is that 
we restricted the dependent variable to a single real line. This indicates to us a significant 
use of restricted complex transformations in converting systems of nonlinear equations into 
their linear analogues. Theorem 3 provides a class of systems of two second-order ODEs 
that can be dealt via complex variables. We now present the basic theorem for complex 
linearization. 

Theorem 3. The necessary and sufficient condition for a system of two second-order ODEs 
of the form 

f" = Ai/'3 _ 2,A^fg'^ _ 2,A^f'^g' + A^g'^ + f'^ - B'g'^ - 2B^f'g' + 

C'f-C\/ + D\ 

g" = SA'f'^g' - A'g'^ + A^f'^ - SA^f'g'^ + 2B'f'g' + B^f' - B^g'^ + 

C^f + C'g' + D^, (38) 

where A^, B\ C*, D\ {i = 1,2) are functions of the variables x, y, f and g, to be solvable by 
complex linearization is that the coefficients satisfy the conditions 

\2Al^ + \2C^A\ - \2AlC'^ - ^A)D^ - ^D^A\ + ^D'^A) - QD^A] + 
l2A^Cl - \2A^Cl + C)^ - Clg + 2C% - 12A^D) - UA^D] + 
12A^D) - 12A^D] + 2B^C} + 2S^Cj - 2B'^C} + 2B'^Cl - SB^Bl + 

SB'Bl-ABlf-ABl^^O, 
\2Al^ + \2C''A\ + \2A?jO^ - ^D'^A) - 6D^Al - 6D^A} + QD^A] + 
\2£Cl + VlA^Cl + C)^ - - 2C}g - 12A^D) - UA^D] - 
12A^D) + l2A^Dl + 2B'^C] + 2B'^Cl + 2B^C} - 2B^Cl - SB'^B^ - 

SB'Bl - AB% + ABlg = 0, 
24:D^Al - 24:D^Al - 6D^B} - ^D^B^ + ^D'^B) - ^D'^B] + 
l2A^Dl - \2A^Dl + ABl^ - AClf - AClg - 6B^D} - GB^D^ + 
GB^Dl - GB^Dl + 3D}f - 3Z5j^ + 6Djg + AC^Cj + AC^C^ - 

AC^Cj + AC^Cl - AC^Bl + AC^Bl = 0, 
2AD^Al + 2AD^Al - GD^Bj - QD^B] - QD^B) + QD^B] + 
l2A^Dl + l2A^Dl + ABl^ - AC^j + AC^g - 6B^D] - QB'^D] - 
GB^Dj + GB^Dl + 3Djf - dD^g - QD]g + AC^C) - AC'^C] + 

AC^C) - AC^C] - AC^Bl - AC^Bl = 0. (39) 
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Proof. Suppose that there exists complex functions, 

A{x, u) = A\x, /, g) + iA\x, f, g), 
Bix,u) = B\xJ,g) + iB\x, f,g), 
C{x,u) = C\xJ,g) + iC^{x, f,g), 

D{x, u) = D\x, /, g) + tD^ix, /, g), (40) 

such that the above system can be mapped into the second-order ODE 

u"{x) = A{x, u)u'''^ + B{x, u)u''^ + C{x, u)u' + D{x, u), (41) 

which is at most cubic in u' and therefore satisfies the necessary condition of linearizability. 
To check the sufficient conditions the set of equations (39) is transformed into the equations 

3A,, + 3A,C + 3Aa, - 3A^D + C^u - 6AD., + BC, - 2BB, - 25,.„ = 0, 

6A,D - 3B^D + 3AD, + B,, - 2C,u - WD^ + 3D„„ + 200^ - CB, = 0, (42) 

which are the Lie conditions. Since the ODE (41) is linearizable, we can obtain its solution 
by the geometric method. This solution can now be written as the pair of real functions, / 
and g. Hence the system of ODEs (38) can be solved by complex linearization. 

To see how the complex variable approach works we present some illustrative examples of 
two dimensional systems of ODEs. 
Examples 

1. Consider the system of ODEs 

//" - 99" = r - 9'^ + {f - 9')wi{x) - 2fgw,{x), 

fg" + gf" = 2/V + 2fgwi{x) + {f - g')w2{x), (43) 

where wi and W2 are arbitrary functions of the time variable x. This was discussed in [12] 
for an anharmonic oscillator with Wi = and W2 = 0. The complex symmetry analysis 
gave interesting insights into the oscillator dynamics. Since the coefficients satisfy the lin- 
earizability conditions therefore the above system is solvable by complex linearization. In 
fact the system is linearizable. The transformation for our purpose is 

^ ^1 T.,/ /■2 , 2N A 1 ........ f 9 



and it reduces system (43) into the linear system of ODEs 

r' = -w,, C=-nJ2. (45) 

X X 

where 

wi=wi{l/x), W2=W2{l/x)- (46) 



2. Two-dimensional Coupled Modified Emden System: 

The modified Emden equation possess numerous dynamical properties in nonlinear oscilla- 
tions. In [8j Chandmsekar et al. explore an important characteristic of such an equation in 
which the frequency of oscillation is independent of amplitude and remains the same as that 
of the linear oscillator. They showed that the amplitude dependence of the frequency is not 
a fundamental property of nonlinear dynamical phenomena. Later, in [TSj, they extended 
the above results and investigated the dynamical properties of N-coupled nonlinear oscillator 
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of Lienard type. Those systems which possess a Hamiltonian structure can be transformed 
into systems of uncoupled harmonic oscillators via contact transformations. We investigate 
the integrability of a system of two cubically semilinear coupled ODEs of Emden type 

/" = -3//' + ^gg' -f + 3fg\ 

g" = -3gf - 3fg' - 3fg + g\ (47) 

The coefficients satisfy the linearizability conditions. Hence the above system is solvable by 
complex linearization. This system corresponds to the cubically semilinear modified Emden 
ODE 

u" + 3uu' + = 0, (48) 

which is linearizable as it satisfies the Lie conditions. This equation arises in many applica- 
tions (see e.g. j24[ [25]). Equation (48) admits the two noncommuting complex symmetries 

d d d 

Zi = |-, Z, = x^-u-^, (49) 
ox ox ou 

which can be used to write down the transformation, 

_ iX-' tXj , V 

X = x--, f/=---, 50 
u 1 u 

that transforms (48) into the free-particle equation \J" = 0. This transformation seems 
odd as X is real while x is complex. The point is that we start and end with a complex 
independent variable restricted to the real line, but in the intervening steps the variable 
moves off it. The procedure is reminiscent of analytic continuation. To check its consistency 
we express the solution of (47) in the new coordinates (x, U) 

U = ax + f3, (51) 

where a and /3 are complex constants. In coordinates {x, u) the above equation yields 

_ 2{x - a) 
x^-2ax-2P' 

which satisfies (48). It generates the solution of the system (47) 

2(x — ai)(x^ — 2aix — 2[5i) + Aa2{oi2X + 132) 



f 



(x2 - 2aix - 2/3i)2 + {2a2X + 2(32^ 



_ 4(x - ai){a2X + (32) - 2a2{x'^ - 2aix - 2/3i) 
^~ (x2 - 2aix - 2/3i)2 + {2a2X + 2(32)^ ' ^ ^ 

which would have been difficult to obtain by other means. Note that here the C-REs are no 
longer preserved by the transformation. 

3. Now consider the Newtonian system of ODEs with velocity dependent forces 

r = 1 + ((/' - x)' - g'^) w, - 2(f - x)g'w2, 

g" = 2{f' - x)g'w, + ((/' - x)' - g") W2, (54) 

where 

wi = wi{2f -x\2g), W2=W2{2f -x\2g). (55) 
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It can be verified that system (54) can be solved by complex linearization as it satisfies the 
criteria of Theorem 3. It corresponds to the complex Newtonian equation 

u" ^1 + {u' - xyw{2u - x'^) (56) 

with quadratic velocity dependent forces. This ODE admits the complex Lie symmetries 

^ d d ^ d r, d 

Zi = l^ + ^7r' 22 = a;— + — . 57 
ox ou ox Ou 

Equation (56) is linearized by the complex transformation 

X = 2u - x^ [/ = X, (58) 

to become 

2U"^-U'w{x)- (59) 
Again (58) is reminiscent of analytic continuation as it is from (real, complex) to (complex, 
real). To check consistency we may take w — 1, i.e., Wi — 1,W2 — 0, to obtain 

2 

M = a + ln2 + y-ln(/3-a:), (60) 

where a and /3 are complex constants. Putting wi = I and W2 = in (54) we get the system 

/" = 1 + + - g'^ - 2xf, 

g" = 2/V - xg', (61) 

with the general solution 

/ = ai-ln2 + ^-lln((/3i-a;)2 + /32), 

g — a2 — arctan ( -r^ — | . (62) 



Also, if we take 

the solution of (56) is 

which gives the solution 



x^ 



f{x) = - + R{x) cos{e{x)), 

g{x) = R{x)sm{9{x)) , (65) 



where 

R{x) 



{ai{pi -x)+ P20i2Y + {(^2012 - a2{Pi - x)y 



2{al + al) 

(32a2 - a2{(3i - x) 

0{x) — 



2{ai{/3i - x) + /32a2y 
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of the system 



, ^ 2^-(f -x)(2/-x2) _ 2^((f -x)2-^-2) 

^ (2/-x2)2+4^2 (2/ _ a;2)2 + 4^2 " ^ ) 

Note that w is an arbitrary complex function that gives a class of systems of ODEs that 
correspond to (56). Thus the linearization of a general equation encodes the linearization of 
a large class of systems of ODEs. 

5. Conclusion and Discussion 

Though the linearizaton procedure for a scalar ODE was fully provided by Lie, there is no 
such complete characterization and procedure available more generally [5l[6l[2T], |24j — [26] ) . 
Geometry gives a procedure that not only provides the invariant characterization but also 
the solution of systems of ODEs [201 EI]- However, it only applies to the class of maximum 
symmetry and it is known that some less symmetric systems of ODEs are also linearizable 
|3S]- A method that retains the power of geometry, but applies to the less symmetric cases 
is needed. 

In this paper we used complex scalar ODEs to write equivalent systems of PDEs and ODEs 
and then required that the original ODEs be linearizable. We provided various examples, 
mainly of mechanical systems, to illustrate the power and use of the method presented. 
For the PDEs we have a guarantee that the C-REs are preserved under the linearizing 
transformation. However, for the ODEs we have no such guarantee. It seems that, when the 
C-REs are preserved, we obtain linearizing transformations and when they are not preserved 
we are not able to linearize the system. However, the solution of the system obtained by 
linearization of the complex scalar ODE is still valid for the equivalent system of PDEs or 
ODEs. 

We called the procedure for ODEs solution by complex linearization. This is an example 
of Penrose's "complex magic" [34]. The solutions are again obtained using the geometric 
linearization procedure. It is hoped that the complex procedures, augmenting the real lin- 
earization for maximally symmetric systems, will give all the classes of linearizable systems 
of ODEs. This line of investigation is being pursued [21 ES] • 

For PDEs we get a systematic procedure to linearize the system and obtain the solu- 
tion. However, it is clear that we do not get all possible solutions. This can be proved by 
considering the linear complex ODE, u" = 0, and writing the corresponding system of PDEs, 

fxx fyy ~l" 2(7x1/ 0, 

gxx - gyy - '^fxy = 0. (67) 

Since there are only two arbitrary (complex) constants that can appear for the scalar ODE 
but infinitely many linearly independent solutions of the system of PDEs, we see that the 
complex linearization procedure cannot exhaust the solutions of the system of PDEs. 

One should be able to extend the complex linearization procedure to third- and fourth- 
order systems by using the results for the corresponding scalar ODEs [ISl 120] straight- 
forwardly. Also the extension to conditional linearizability of systems could be obtained 
[27 ] [25 ] E5]. However, in the former case the power of the geometric approach would be lost. 
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In the latter case, if the root equation is second-order, the geometric method would provide 
the solution of the system. 

A bigger problem is the extension to higher- dimensional systems. The extension to di- 
mensions of 2n can be obtained by iterative use of the complex method. Starting with a real 
system of n-dimensions and converting to complex variables, we could get 2n-dimensional 
systems as well. However, it would appear that the iterative procedure and the complexifi- 
cation of the system would yield different linearizable classes in general. Furthermore, this 
would not provide a means of deahng with systems of odd dimensions. It would be inter- 
esting to explore the various ramifications of extension of complex linearization to higher 
dimensions. 
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